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Abstract 

We prove that the purely imaginary square well generates an infinite number of bound 
states with real energies. In the strong-coupling limit, our exact VT symmetric 
solutions coincide, utterly unexpectedly, with their textbook, well known Hermitian 
predecessors. 
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1 Introduction 



An interest in the imaginary cubic anharmonic oscillators dates back to their per- 
turbation analysis by Caliceti et al The simplified homework example with the 
mere two-term non-Hermitian Hamiltonian 

Hbz =p'^ + ix^ 

has been proposed by D. Bessis and J. Zinn-Justin who had in mind its possible 
applicability in the context of statistical physics [Q. The example has been re- 
vitalized by C. Bender et al due to its further possible methodical relevance in the 
relativistic quantum field theory 0]. 

The apparent reality of the spectrum of energies Ebz proved quite puzzling and 
inspired a conjecture of the existence of the whole "modified quantum mechanics" 
paying attention to similar Hamiltonians [^]. This conjecture opened many new and 
interesting questions when it replaced the current Hermiticity of Hamiltonians by 
a weaker condition of their commutativity with a product VT of the spatial parity 
V and the complex conjugation T where the latter factor is to be understood as a 
one-dimensional version of the operator of time reversal. 

The recent discussion between C. Bender and A. Mezincescu pointed out that 
one of the key problems of the new studies lies in the ambiguity of the spectrum 
which depends quite crucially on our choice of the boundary conditions which can 
be, in general, complexified 0. The fragile character of the reality of energies has 
been confirmed by the WKB and perturbative studies and by the quasi-exact and 
exact models ||^ where the admissible unavoided level crossings ^ prove sometimes 
followed by the spontaneous breakdown of the VT symmetry |[T0[| . 

In such a context we intend to propose an extremely elementary VT symmetric 
model which would replace the BZ interaction ix^ (which admits just a numerical 
treatment) by its exactly solvable square-well analogue. 
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2 Model 



In a search for analogies between the solvable and unsolvable models in one dimen- 
sion, all the possible forms of a confining well are often being approximated by the 
ordinary real and symmetric square well 

X G (-oo,-7r)U(7r, oo), 
0, X G (— vr, tt). 

In this spirit we intend to replace here the above-mentioned antisymmetric and imag- 
inary homework potential Vbz{x) = ix^ hj its elementary square- well analogue 

—iT"^, a: G (— oo, — tt), 
0, X E (— TT, 7r), 

X G (tt, oo). 

Schrodinger equation which appears in such a setting, 



2m dx^ 



ip{x) = Eip{x) 

will be complemented here by the standard L^(R) boundary conditions 



(2) 



?/'(±oo) 



0. 



(3) 



The well known VT symmetric normalization convention will be employed, with a 



free real parameter G in the unbroken PT-symmetry requirements |11 



7/^(0) = 1, d,,ij{^)=iG. 



(4) 



Putting h = 2m = 1 and using the ansatz 



(5) 



coskx + Bsinkx, xG(0,7r), k'^ = E, 

(L + iN) exp{-ax), xG(7r, oo), a"^ = i T"^ - k"^ , 

we guarantee its full compatibility with the symmetry requirements by the choice 
of the purely imaginary constant B = iG/k'm the wave functions (0). 
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3 Matching conditions at x = n 



We may split a = p + iq in its real and imaginary part with p,q > 0. This gives the 
rules p"^ + k"^ = and 2pq = as a consequence. They are easily re-parameterized 
in terms of a single variable a, 

T 



p = q cos a, k — q sina, q — 



a e (0,7r/2). 



(6) 



-\/2cosa' 

In this language the standard matching at the point of discontinuity is immediate, 
cosk7r + B sinkn — (L + iN) cxp(— cttt), 
— sinkn + B coskn — —^(L + iN) exp(— cttt). 

fV 

After we abbreviate a/k — — tan fin, we get an elementary complex condition of the 
matching of logarithmic derivatives at x — tt, 



G — —ik tan{k + fl)'!r. 



(7) 



Its real part defines our first unknown parameter, G = G{a). Due to our normaliza- 
tion conventions, the imaginary part of the right-hand-side expression must vanish, 
Re[tan(/c-|-Jl)7r] = 0. An elementary re-arrangement of such an equation acquires the 
form of an elementary quadratic algebraic equation for X — tan kir. Its two explicit 
solutions read 



^1 = 



p+_q 
k ' 



P 



k 



or, after all the insertions. 



tan 



tan 



vrTsino;*^"'"^ 
V2 cosa(+) 

ttT sino;'^") 



tan 



TT — a 



(+) 



V2 



cos a 



(-) 



tan 



(-) 



(8) 



(9) 



(10) 



These equations specify, in implicit manner, the two respective infinite series of the 
appropriately bounded real roots a — a!;^'' e (0,7r/2). 
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4 Energies 

Even before any numerical considerations we immediately see that for a G (0,7r/2) 
the left-hand-side arguments [. . .] in eqs. (P) and (|10D run from zero to infinity. 
Their tangens functions oscillate infinitely many times from minus infinity to plus 
infinity. Within the same interval, the limited variation of the argument a makes 
both the eligible right-hand side functions monotonic, very smooth and bounded, 
tan[(7r — a;*^+^)/2] G (1, oo) and tan[a;'^~Y2] G (0, 1). A priori this indicates that our 
roots k = k{al^^) will all lie within the fairly well determined intervals, 

ki+^ e (^n+^,n+^y n = 0,l,..., 

k^^^ G (^m + ^, m + 1^ m = 0, 1, 

After such an approximate localization of the roots, an unexpected additional merit 
of our parametrization manifests itself in an unambiguous removal of the tangens 
operators from both eqs. @ and ([l^)- This gives the following two relations, 

ki^'=n+'--^. e^=m + l-i^, a;i^) = ^ G (0,1). 

/ 4 4 TT 

After an elementary change of the notation with cu!^^ = uj2n and tu!^^ = co2n+i, we 
may finally combine the latter two rules in the single secular equation 



sin (l^jv) = " Y ^ (i^^) N = 0,1,..., (11) 

In a graphical interpretation this equation represents again an intersection of a 
tangens-like curve with the infinite family of parallel lines. This is illustrated in 
Figure 1. The equation generates, therefore, an infinite number of the real roots 
ujy G (0, 1) at all the non-negative integers = 0, 1, . . .. 
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5 Wave functions in the weak coupling regime 

Equation (|^ in combination with eqs. (P) and (p!oD determines the real parameter 

G = = — (12) 

q ±p 

responsible for the behaviour of the wave functions near the origin [remember that 
B = iG/k'ui eq. (|^)]. For its deeper analysis let us first introduce an auxiliary linear 
function of uj and A^, 



2N + 2~ujn ^ (N+l/2 N+V 



AT V 2T ' 2T ; 

and re-interpret our secular eq. (|lTD as an algebraic quadratic equation with the 
unique positive solution, 

cos {'^^^n] = 1 — ■ (13) 

This is an amended implicit definition of the sequence un- As long as the right hand 
side expression is very smooth and never exceeds one, the latter formula re-verifies 
that the root cua? is always real and bounded as required. 

In the domain of the large and almost constant 1 (i.e., for the small square- 
well height T or at the higher excitations), our new secular equation (^) gives a 
better picture of our bound-state parameters un = i — Vn which all lie very close to 
one. The estimate 

vr 115 

— 'nN = arcsin , ^ + . . . 

2 R + y/WTT 2R 48i?3 ^ 

represents also a quickly convergent iterative algorithm for the efficient numerical 
evaluation of the roots ujn. One can conclude that in a way compatible with our a 
priori expectations, the value of p = pn = Reo" ~ q/2R is very close to zero and, as a 
consequence, the asymptotic decrease of our wave functions remains slow. We have 
q = Qn = Inio" ^ k so that, asymptotically, our wave functions very much resemble 
free waves exp{—ikx) . In the light of eq. (|1^) we have also tp{x) ~ exp{—ikx) near 
the origin. 
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6 Wave functions in the strong coupling regime 

For the models with a very small R (i.e., for the low- lying excitations in a deep well 
with T ^ 1) we get an alternative estimate 



TT 1 

- cjat = arcsm W - 



R - i Vl + _ 1 



In the limit i? ^ the present spectrum of energies moves towards (and precisely 
coincides with) the well known levels of the infinitely deep Hermitian square well of 
the same width I = {—tt,tt) (cf. eq. (|I]) with 5* oo). In this sense, the "complex- 
rotation" transition from the Hermitian well V^^^^\x) of eq. (|1]) (with S ^ 1) to 
its present non-Hermitian VT symmetric alternative V^^^'^^\x) of eq. (|]) (with 
T ^ 1) proves amazingly smooth. 

The wave functions exhibit the similar tendency. In the outer region, they are 
proportional to exp(— pa;) and decay very quickly since p = 0{R~^^'^). The parameter 
G^^^ becomes strongly superscript-dependent, 

G(+) = — ^ = 0{R'/^), G(-) = -{q + p) = 0{R-'/''). 
q + p 

This means that in the interior domain of a; G (— vr, it), the wave functions with the 
superscript and become dominated by their spatially even and odd components 
cos kx and sin kx, respectively. In this sense, the superscript mimics (or at least keeps 
the trace of) the quantum number of the slightly broken spatial parity V. 

We can summarize that our present VT symmetric model is, unexpectedly, quite 
robust. Almost irrespectively of the coupling T, the spectrum is unbounded from 
above and remains constrained by the inequalities 

(7V + V2)^ < < (E±ll. (M) 

The analogy between our exactly solvable square-well model and the standard or 
"paradigmatic" VT symmetric Hamiltonian Hbz appears to be closer than expected. 



6 



7 Outlook 



The exact solvability of our present purely imaginary square well model throws a 
new light on some properties of the VT symmetric wave functions which are hardly 
accessible by approximative techniques. In the nearest future, one can expect that the 
further detailed study of the "PT symmetric square wells will give new answers to the 
recent puzzles as formulated in ref. [|1^] and concerning the irregular behaviour of the 
nodal zeros in the complex plane. Our present example indicates that a surprising 
alternative to the Sturm Liouville oscillation theorem could, perhaps, emerge in 
connection with the study of zeros of the separate real and imaginary parts of the 
VT symmetric wave functions. 
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Figure captions 

Figure 1. Graphical solution of eq. ([TT| ) (?/ = un/'^-, T = 1) 
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Figure 1 




